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ABSTRACT 

Aims. RXS J113155.4-123155 (z=0.66) is a quadruply imaged lensed quasar with a resolved Einstein Ring. The goal of this paper is to provide 
a full characterization of this system, and more particularly accurate astrometry and photometry. These observational constraints constitute 
a mandatory ingredient for the precise determination of the lens mass profile, the derivation of the Hubble constant Ho from time delay 
measurements and investigations on the presence of massive substructures in the lensing galaxy. 

Methods. Visible and near-infrared imaging observations of RXS J113155.4-123155 were carried out at various epochs using several ground 
based telescopes and the HST The frames have been deconvolved using the MCS algorithm. A Singular Isothermal Ellipsoid (SIE) + external 
shear has been used to model the lensing galaxy potential. 

Results. MCS deconvolution enables us to separate the flux of the QSO (point-like images) from that of its host galaxy and to accurately 
track the flux variations of the point-like images in various filters. The deconvolved frames unveil several multiply imaged structures in the 
Einstein ring and an unidentified object in the vicinity of the lensing galaxy. We discuss the lightcurves and the chromatic flux ratio variations 
and deduce that both intrinsic variability and microlensing took place during a span longer than one year. 

We demonstrate that microlensing may easily account for the so called anomalous flux ratios presented in the discovery paper. However, the 
observed flux ratios are still poorly reproduced when modeling the lens potential with a SlE+shear. We argue that this disagreement can hardly 
be explained by milli-lensing caused by substructures in the lensing galaxy. A solution proposed in Paper 11 consists in a more complex lens 
model including an octupole term to the lens gravitational potential. 
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1. Introduction 

RXS J113155.4-123155 (hereafter J1131) is one of the nearest 
confirmed multiply imaged AGN. The source at Zs = 0.658 is 
lensed by an elliptical galaxy at zi — 0.295 (Sluse et al. 120031 
hereafter Paper T). This system is a long axis quad with an im¬ 
age configuration very similar to B1422-1-231 (Patnaik et al. 
fT^ : three merging images (B-A-C; typical of a source ly¬ 
ing close to a cusp caustic) face the faint saddle-point image 
D lying close to the lensing galaxy (G). For such a system (i.e. 
cusp configuration system), the magnification behaviour is well 
understood - the flux of the middle image should be equal to 
the total fluxes of the two outer images (Schneider & Weiss 
119^ . However, this generic prediction was strongly violated 
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at the epoch of the discovery, suggesting the likely presence 
of small-scale structures in the lensing galaxy (Keeton, Gaudi 
& Fetters IMBl l2005l hereafter KGP03, KGP05). Thanks to 
the data set presented here, we re-examine the values of the 
flux-ratios with regard to their temporal and chromatic varia¬ 
tions. Additionally, because this quad is quite bright and shows 
lensed images with a wide angular separation {29e ~ 3.6"), it 
is potentially a target of interest for time delay measurements. 
Also, the bright Einstein ring detected from the optical to the 
near-infrared offers unique constraints on the lens modeling 
(Kochanek et al. l2001l l that may break the degeneracy between 
the intrinsic lens galaxy shape and the external shear (Keeton 
et al. I2000t . However, the bright Einstein ring introduces an 
additional complexity to the photometry measurement because 
it is superimposed over the lensed point-like images. This mo¬ 
tivates our extensive discussion of the sources of systematic 
photometric errors affecting the individual lensed QSO images 
and of the photometric accuracy that can be reached for J1131. 

The observations presented here constitute the first part of 
a follow-up study of this system that also includes detailed lens 
models and source reconstruction from HST imaging. X-ray 
imaging as well as optical and NIR spectroscopy. Sects. E&El 
describe the data reduction and analysis. Sectionincludes a 
description of the ring morphology, a simple modeling of the 
lensing galaxy and a discussion of the flux (time and chromatic) 
variability of the lensed images. Einally, Sect.|3summarizes the 
main results and presents our conclusions. 

We adopt throughout this paper Hq - 65 km s ' Mpc ', 
Qo = 0.3 and Aq = 0.7; magnitudes are computed in the Vega 
system. 

2. Observations and reductions 

In this section, we present NIR and optical imaging ob¬ 
tained for J1131. Near Infrared observations include images 
recorded with ESO telescopes, with CFHT and HST. The im¬ 
ages recorded at visible wavelengths have been obtained with 
the FORSl and FORS2 instruments at the Cassegrain focus of 
Antu (UTl) and Yepun (UT4) in Paranal. Tabled summarizes 
the technical data of the various instrumentations used and lists 
the observational characteristics of each data set. 

2.1. Ground-based NIR observations 

First, J band images of J1131 have been obtained with SOFI. 
The conditions were photometric and 6 standard stars have 
been observed during the night enabling us to derive the magni¬ 
tude of an object producing 1 count/s of 23.15 + 0.01 (hereafter 
zero-point). 

Secondly, images of J1131 have been obtained in the Ks- 
band with the Short Wavelength Imaging camera of the ISAAC 
instrument. Thanks to the short DIT (10s) used, the combined 
flux from the sky and QSO images in any pixel never ex¬ 
ceeded 10000 ADU, a level below the non-linearity regime of 
the ISAAC detector. Conditions were photometric and a zero- 
point of 24.17+0.07 has been deduced from the standard stars 
observed during the night. 


Finally, a set of NIR images of J1131 has been obtained at 
a single epoch with CFHT-IR through the broad band 7#5133, 
//#5209 and K'#5337 biters. 

For the whole set of data, standard NIR reduction proce¬ 
dures were applied to subtract the dark and to bat beld the 
images using a normalized bat beld. Sky subtraction and co¬ 
addition of the reduced dithered frames were performed us¬ 
ing the xdimsuin IRAF* package. The residual sky pattern has 
been btted and subtracted using the SExtractor software v 2.3.2 
(Bertin,IT996j. 

2.2. Ground-based optical observations 

Observations with the FORS2 instrument have been obtained 
with the High Resolution (HR) collimator and the CCD detec¬ 
tor in the 2x2 binning mode. Conditions were not photometric. 

B, V, R, I images of J1131 were obtained in April 2004 
under photometric conditions with the FORS1 instrument and a 
CCD in the unbinned mode. Zero-points for the Bessel B, V, R, 
I biters were found to be 27.76+0.02, 28.12+0.02, 28.06+0.04 
and 27.25+0.01 respectively. 

Standard reduction procedures were applied to the data in¬ 
cluding bias subtraction and bat belding with a normalized bat. 
The photometric measurements of the standard stars have been 
performed with the phot IRAF task. The sky has been sub¬ 
tracted using the SExtractor software v 2.3.2 (Bertin. n996> . 

2.3. HST FI 60W observations 

HST F160W observations of J1131 have been obtained on 
November 17, 2003 with the Hubble Space Telescope (HST- 
GO-9744; PI; C.S. Kochanek) as part of the CfA-Arizona 
Space Telescope Lens Survey (CASTLES). Out of the eight 
dithered frames of J1131 obtained with the NIC2, two frames 
have been excluded from our analysis because of an error bag 
associated with pixels located in the brightest point-like im¬ 
ages. We reduced the images using the IRAF calnica and 
calnicb tasks as recommended in the NICMOS data hand¬ 
book (Dickinson et al. I2002t . The bux calibration (Dickinson 
et al. 1200^ has been derived from PHOTNU=1.4981 10 ^ 
Jy sec DN * and from the averaged bux density of VEGA 
< Fy >- 1043.5 Jy (June 2004 values). 

3. Data analysis 

3.1. The method 

The Einstein ring of J1131 is the brightest one known to¬ 
day. Because it is superimposed over the lensed QSO images, 
standard photometric methods like PSF btting may bias the 
point-like bux measurements. We thus measured relative pho¬ 
tometry using the MCS deconvolution method (Magain et al. 
which proved to be very efficient in measuring buxes 

’ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Log of the observations. Epoch = epoch number as used in the text. CONAD = electron to ADU conversion factor. 
NExp = number of dithered frames. Exp = total integration time of a single frame. Eor the NIR data, Exp = NDITxDIT, where 
NDIT is the number of sub-frames with an integration time of DIT seconds. Only the total exposure time per frame is reported 
for the NIC2 images. EWHM = mean seeing measured on frames. Note that the SOEI pixel size has been re-measured to be 
0.289", superseding the published ESO value of 0.292". 


Date 

(dd-mm-yyyy) 

Epoch 

Instrument 

(Telescope) 

RON 

(e^/pixel) 

CONAD 

(e-/ADU) 

Fix. size 

(") 

Filter 

NExp 

Exp 

(s) 

FWHM 

(") 

21-11-2002 

1 

SOFI (NTT) 

11 

5.3 

0.289 

J 

24 

2x30 

0.7-0.9 

12-04-2003 

2 

ISAAC (UTl) 

11 

4.6 

0.1484 

Ks 

48 

6x10 

0.3-0.55 

21-04-2003 

3 

FORS2 (UT4) 

4.2 

1.25 

0.1263 

R special 

4 

80 

0.65-0.85 

02-05-2003 

4 

FORS2 (UT4) 

4.2 

1.25 

0.1263 

R special 

6 

40 

0.8 

26-05-2003 

5 

FORS2 (UT4) 

4.2 

1.25 

0.1263 

B Bessel 

8 

497 

0.65-0.9 

26-05-2003 

5 

FORS2 (UT4) 

4.2 

1.25 

0.1263 

B Bessel 

7 

125 

0.75-1.0 

17-06-2003 

6 

FORS2 (UT4) 

4.2 

1.25 

0.1263 

R special 

6 

40 

0.45-0.50 

17-06-2003 

6 

FORS2 (UT4) 

4.2 

1.25 

0.1263 

V Bessel 

20 

59 

0.5-0.65 

17-11-2003 

7 

NIC2 (HST) 

26 

5.4 

0.075 

Fmw 

5 

640 

0.12 

17-11-2003 

7 

NIC2 (HST) 

26 

5.4 

0.075 

Fmw 

3 

704 

0.12 

18-11-2003 

7 

FORS2 (UT4) 

4.2 

1.25 

0.1263 

B Bessel 

3 

240 

1.0 

18-11-2003 

7 

FORS2 (UT4) 

4.2 

1.25 

0.1263 

V Bessel 

4 

120 

0.8-1.0 

18-11-2003 

7 

FORS2 (UT4) 

4.2 

1.25 

0.1263 

R special 

4 

80 

0.85-0.95 

09-02-2004 

8 

CFHT-IR (CFHT) 

15 

2.35 

0.211 

y#5133 

4 

6x120 

0.6-0.8 

09-02-2004 

8 

CFHT-IR (CFHT) 

15 

2.35 

0.211 

H#5209 

4 

6x60 

0.6-0.8 

09-02-2004 

8 

CFHT-IR (CFHT) 

15 

2.35 

0.211 

K'#5337 

8 

3x60 

0.55-0.9 

12-04-2004 

9 

FORSl (UTl) 

5.2 

1.61 

0.100 

B Bessel 

3 

60 

0.7-0.8 

12-04-2004 

9 

FORSl (UTl) 

5.2 

1.61 

0.100 

V Bessel 

4 

30 

0.65-0.85 

12-04-2004 

9 

FORSl (UTl) 

5.2 

1.61 

0.100 

R Bessel 

4 

20 

0.7-0.85 

12-04-2004 

9 

FORSl (UTl) 

5.2 

1.61 

0.100 

I Bessel 

4 

20 

0.6-0.7 


of point-like components superimposed over a varying back¬ 
ground (Burud et al. 119981 . Thanks to this code, images are 
deconvolved to an improved but finite resolution using a ker¬ 
nel constructed from PSE stars (Courbin et al. 119981 . Eor HST 
data, the kernel is constructed based on the PSE generated with 
the Tiny Tim v 6.1 software (Krist & Hook l200^ . In the other 
cases, the kernel is constructed from relatively isolated bright 
held stars. Additionally, since the co-addition of different im¬ 
ages often degrades the PSE, we generally worked on the indi¬ 
vidual frames obtained at different dithered positions. 

Instead of deconvolving each image separately, we used a 
modihed version of the code that simultaneously deconvolves 
a set of individual frames taken with the same instrumental 
setup (Burud et al. 120001 . leading to better constraints on the 
astrometry of the quasar images (an output of the deconvo¬ 
lution process), on the shape of the Einstein ring and of the 
lensing galaxy, using the S/N of the whole data set. Eor each 
data set, deconvolved frames are reconstructed to a common 
improved resolution. Technically, the relative positions of the 
point-like components are htted on the reference frame while 
their flux and the frame offset with respect to the reference vary 
independently for each frame. Eor the data sets with the largest 
pixel size (i.e. CEHT-IR and SOEI), MCS could not fit accu¬ 
rately both the offset with respect to the reference frame and 
the point-like image positions. In that case, the latter have been 
fixed to the ones derived from the HST observations. 

In the following, we have simultaneously deconvolved all 
the frames; i) for each epoch (in a given hlter) separately; ii) for 
all the epochs available (in a filter). Single epoch deconvolution 


is best suited to retrieve absolute photometry (Sect. l3..3~n since 
it is not affected by the flux scaling between several epochs. It 
is also the only method applicable to the I band and NIR bi¬ 
ters (i.e. one epoch data). Multi-epoch deconvolution is more 
appropriate to retrieve accurate bux ratios (Sect. and 

to measure the temporal bux variations (Sects. msi&Eis 
Typical photometric random errors obtained with MCS reach 
the photon noise uncertainty. 

Although the MCS deconvolution enables one to retrieve 
accurate photometry, the results are affected by systematic er¬ 
rors (Burud et al. mainly associated with the regulariza¬ 

tion term A. The latter works like a local smoothing term by re¬ 
ducing the high frequencies of the background. Thus, depend¬ 
ing on the choice^ of A, a variable fraction of the background 
bux is included in the point-like bux measurements. We se¬ 
lected realistic values of A based on the batness of the resid¬ 
ual map (i.e. the difference between the real image and the de¬ 
convolved image re-convolved with the kernel, and divided by 
the standard deviation in each pixel of the real image). More 
specibcally, the higher level of regularization is picked when 
we observe more than 3 cr standard deviation in the residual 
map (regularization is swamping -true- structures existing in 
the original image) and the lower level is chosen when impor¬ 
tant Gibbs oscillations appear in the deconvolved frame. In the 
following, we run MCS with three acceptable levels of regular¬ 
ization; a lower and upper bound -resp. Tiow and dhigh- and an 
intermediate value dbest considered as the best value. 

^ A high value of A means a small regularization and a low value 
means a strong smoothing of the deconvolved image. 
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Table 2. Relative positions of the different lensed components 
(BCD), the lensing galaxy (G) and the companion object (X) 
with respect to A. These are deduced from a gaussian fit¬ 
ting (except X with MCS) on the HST NIC2 image (F160W) 
and from the MCS deconvolution of the FORS2 and ISAAC 
frames. The last column lists the Icr random error in both right 
ascension and declination. For the ground-based data, the error 
on the relative position is reported. Icr systematic errors related 
to scaling, orientation and distortion uncertainties are estimated 
to be 0.003” (see text for details). 


ID 

Aa cos5(") 

Ah(") 

lo-n 


F160W 

ground 

F160W 

ground 

F160W 

ground 

A 

0 

0 

0 

0 

0.001 

- 

B 

0.032 

0.030 

1.184 

1.190 

0.001 

0.001 

C 

-0.590 

-0.589 

-1.117 

-1.117 

0.001 

0.002 

D 

-3.115 

-3.120 

0.875 

0.881 

0.001 

0.002 

G 

-2.027 

-2.027 

0.607 

0.604 

0.001 

0.002 

X 

-1.936 

- 

1.100 

- 

0.011 

- 


Another source of systematic photometric errors may come 
from an imperfect PSF model. However, this appears as devia¬ 
tions from a null residual map and the associated photometric 
bias is corrected by means of the comparison of deconvolu¬ 
tion and aperture photometry of held stars (see e.g. Sect. 13.3.11 . 
Note that in such cases, since the residual map cannot be used 
to pick up /flow and dhigh, we select them based on the visual in¬ 
spection of the deconvolved frame: for a too small value of A, 
the background does not look like a ring-i-lens galaxy, but just 
like a very diffuse background while we determine dhigh via the 
appearance of important Gibbs oscillations. 

By default, the error bars reported in this paper are Icr stan¬ 
dard errors on the mean. 

3.2. Astrometry 

The relative astrometry is derived from the HST data using two 
methods. First, we performed 2D gaussian htting of the individ¬ 
ual lensed images and of the lensing galaxy (each component 
isolated in a box of 13x13 pixels) with the i2gaussfit IRAF 
task. Mean positions and associated errors are calculated from 
the measurements on the 6 individual frames. 

In the second method, the relative astrometry of images B, 
C, D and of the lens galactic nucleus G relative to A is derived 
simultaneously using MCS. Unfortunately, this method suffers 
from the poor PSF sampling and is very sensitive to the prior 
knowledge of the relative offset between the frames. Hence, 
we have used the gaussian htting results to hx the mean rel¬ 
ative offset and measured the mean position of object X (see 
Sect. l4.lt with MCS. Note that MCS positions of A, B, C, D 
and G are in statistical agreement with the gaussian htting re¬ 
sults. 

The measured positions have been converted to the 
RA/DEC system using a plate scale^ of 0.07588"/pixel in x 


and of 0.07537 "/pixel in y and a rotation angle of 59.202° (E 
ofN). 

Relative astrometry is also derived from the best ground- 
based observations (i.e. EORS2 and ISAAC data). Since the 
spatial sampling of the ground-based PSEs fully satishes the 
sampling theorem, MCS leads to an accurate relative astrom¬ 
etry. The relative positions of the lensed quasar images and of 
the lensing galaxy found in each band (namely B, V, R and Ks) 
do agree between each other, within typically less than 0.01". 
The standard deviation of the positions measured in each hlter 
is used to derive the Icr standard error on the mean relative po¬ 
sitions. However, because of the small separation of the triplet 
of images A, B and C, the measured positions of B and C are 
likely correlated with that of A, resulting in an underestimate 
of the error. Consequently, we have arbitrarily and quadrati- 
cally added a 0.001” error to the Icr standard errors of B and 
C. 

The systematic error on the position is estimated to be 
0.003". This error is typical of the systematic error that may 
affect the relative astrometry of small separation systems on 
NIC2 data (Impey et al.[T998). Eor the ground-based data, the 
systematic error is inferred from the recalibration of the pixel 
size of the EORS frames with the MIDAS implementation of 
posl (Walter & West 1 198^ and with GAIA^ (v 2.5.3). Table|2 
summarizes the relative positions we derived for J1131 based 
on the HST and ground-based data. Both sets of results do 
agree within 0.006". 

3.3. Photometry 

3.3.1. Absolute photometry 

We derive absolute photometric measurements with MCS in 
the Bessel B, V, R, I bands and for the J, F\60W and Ks fil¬ 
ters (Table |3 based on the fluxes obtained with Tbest- Sets of 
individual frames are used except for the 7-SOEI (Ks) filters 
where a set of 4 (8) stacked frames of 6x1 minutes is con¬ 
sidered. The peak to peak systematic variations introduced by 
A are estimated from the results obtained for diow and Thigh. 
Eluxes are converted into the Vega magnitude system using 
the zero-points reported in Sect. 12. H and applying an extinction 
correction (see the ESO webpages for extinction coefficients). 
Additionally, for each set of ground-based images, we com¬ 
pared the aperture photometry obtained with the phot IRAE 
routine and SExtractor to MCS photometry of at least 4 stars in 
the field. All the measurements are consistent, suggesting the 
absence of systematic errors caused by the estimate of the PSF 

For the HST NIC2 data, the Tiny Tim PSF poorly re¬ 
produces the observations. Consequently, we find several ac¬ 
ceptable solutions for each possible value of A (Thest but also 
Tiow.high) suggesting that the minimization routine in MCS is 
trapped in one of the numerous local minima. Hence, the ran¬ 
dom error bars of the F160W photometry presented in Tabled 
are likely underestimated and the associated peak to peak vari- 


* GAIA is a derivative of the Skycat catalog and image display tool, 
developed as part of the VLT project at ESO. Skycat and GAIA are 
http: //WWW . stsci . edu/hst/nicmos/performance/platescafee softwares under the terms of the GNU copyright. 
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Table 3. Photometry in the B-V-R-I Bessel filters (epoch 9), in the J band (epoch 1), in the //-band {F\60W Vega normalized; 
epoch?) and in the /Tsband (epoch?) for images A, B, C and D of J1131 as obtained with MCS. The notation is t«best±o‘randoin -sys’ 
where /Wbest is the mean value obtained for dbesu CTrandom is the quadratic sum of the 1 cr error on the mean, the error on the extinction 
coefficient and the error on the zero-point (this one amounts to 0.03 mag for NICMOS photometry) and +sys is the peak to peak 
variation between extreme values of A. Last column gives the FWHM (in arcsec) of the PSF in the deconvolved frame. When 
null error bars appear in the table, the error bar is of only a few milli-magnitudes. 


Filter 

epoch 

A 

B 

C 

D 

FWHM 

B 

V 

R 

I 

J 

/’1601V 

Ks 

9 

9 

9 

9 

1 

7 

2 

18.09 ± 0.02/"™ 
17.88 ± 0.02/°;“ 
17.76 ± 0.04/°;“ 
17.43 ± 0.01/°;“ 
16.83 ± 0.03/°;“ 
15.80 ± 0.04/°;“^ 
14.72 ± 0.07/°;“* 

18.14 ±0.02/";°° 
17.83 ± 0.02/°;“ 
17.66 ± 0.04/°;“ 
17.42 ± 0.01/°;“ 
16.49 ± 0.03/°;“ 
16.08 ± 0.03/°;“ 
15.11 ±0.07/°;°° 

19.19 ±0.02/°;“ 
18.96 ± 0.02/°;“ 
18.77 ± 0.04/°;J° 
18.44 ± 0.01/°;J^ 
17.48 ± 0.04/°;“ 
16.95 ± 0.04/°;2‘ 
15.72 ± 0.07/°;“ 

20.64 ± 0.02/";“ 
20.41 ±0.02/°;* * 
20.21 ± 0.04/°;“ 
19.72 ± 0.01/°;“ 
19.00 ± 0.04/°;“ 
18.36 ± 0.07/°;“ 
17.17 ±0.07/°;“ 

0.30 

0.30 

0.35 

0.30 

0.29 

0.08 

0.15 


ations may also be biased. Although we do not know by which 
amount, we trust these results since they are found in good 
agreement with independent PSF fitting photometry presented 
in Claeskens et al. (accepted; hereafter Paper//). 

3.3.2. Relative photometry 

We report in this section the flux ratios between the lensed im¬ 
ages deduced from the multi-epoch deconvolution when possi¬ 
ble and from single epoch data otherwise (cf Sect. 13.11 . 

Multi-epoch deconvolution consists in deconvolving simul¬ 
taneously a set of frames obtained in a Alter at several epochs. 
We report the results obtained for B-, V- and P-band data^ at 
3, 3 and 5 epochs. For all epochs and all bands, except epoch 6 
in the V band, all individual frames have been used. For epoch 
6 in the V band, we only included the 5 best-seeing frames. For 
each Alter, we rebinned the FORS1 frames to the sampling of 
the FORS2 frames. Then, we applied frame to frame zero-point 
correction calculated based on photometric measurements of 
common stars in each held. Then, we ran MCS on the regis¬ 
tered data and derived the magnitude differences and flux time 
variations. Due to a small apparent systematic bias (~0.01 mag) 
of the MCS photometry, we corrected our measurements based 
on the comparison of the aperture photometry and MCS pho¬ 
tometry of 4 held stars. Final results are presented in Table 0] 
(11 first rows). 

The reduction procedure of single epoch data has already 
been described in Sect. 13. 3. f1 Only CFHT-IR data have not yet 
been discussed. For those, the complete set of individual frames 
have been simultaneously deconvolved with a Anal FWHM of 
0 . 21 ". 

The full results are reported in Table 0] The photometric 
(random) errors were calculated based on the dispersion of the 
measurements obtained for the different frames. Following this 
method, the Poisson noise and the error coming from the im¬ 
perfection of the kernel used for deconvolution are included in 

^ Although J (K) hand data exist at two epochs, the filter response 
curves are significantly different. Flux variations are thus less reliahle 
and not reported. 


the error bars. This is confirmed by the fact that estimated er¬ 
rors are larger or equal to photon noise. The systematic error 
was estimated from the results obtained for diow and dhigh (cf 
Sect. 1.3. H . We observe, between these two limits, typical varia¬ 
tions of the flux ratio by 5 up to 15% from the V to the K bands 
while it is < 2% in the B band for which the contrast between 
the ring and the lensed QSO images is maximum. At a given 
wavelength, these systematic errors are reduced when the S/N 
in the Einstein ring increases. 

Finally, note that whenever possible, we have compared the 
flux ratios deduced during the multi-epoch process to single 
epoch deconvolution results. Both measurements are basically 
identical® but multi-epoch results are less sensitive to system¬ 
atic errors because of the increased S/N in the ring that im¬ 
proves the separation between point-like and extended flux. 

3.3.3. Caveat 

Because of the presence of the Einstein ring superimposed over 
the quasar images, one may ask whether the measured point¬ 
like fluxes are associated with the quasar only (i.e. compact 
emission) or whether they are significantly mixed with emis¬ 
sion originating from a more spatially extended (partly unre¬ 
solved) region of the host galaxy. In the following we examine 
how the core of the host galaxy may modify the observed fluxes 
of images A, B and C. 

The observed flux ratios between images A, B, C are con¬ 
taminated in two ways by the ring. 

The first one is a purely instrumental effect. The resolved 
ring plays the role of a spatially variable background, whose 
intensity is increasing below the unresolved PSF because this 
region corresponds to the innermost parts of the host. However 
the deconvolution process tends to produce a ’’flat” background 
below the PSF whose intensity closely matches the value in 
its close vicinity. Since the surface brightness of the resolved 
background ring is the same around each point-like image, the 
flux of the faint images (i.e. weakly amplified) is more contami- 

® Optical single epoch magnitude difference can be inferred for 
epoch 9 in Tabled 
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Table 4. Magnitude differences (Sect. 13.3~^ . The notation is mtest ± cTrandom -syL where tWbest is the mean value obtained for dbest^ 
CTrandom is the 1CT eiTor on the mean and +sys refers to the peak to peak variations between extreme values of A. Col 1: filter; Col 
2: epoch for which the magnitude differences are given; Col 3: R/j(A) = /,7/y for image A. i is the reference epoch (bold face) for 
the filter considered while j is the epoch given in Col. 2. Col. 4-5-6: average magnitude differences between image A and resp. 
B, C, D. Col. 7-8: idem between image B and resp. C, D. Col. 9: FWHM (in arcsec) of the PSF in the deconvolved frame. Since 
the flux variations (Col. 3) are not affected by the regularization, systematic error bars are not expressed. When null error bars 
appear in the table, the error bar is only of a few milli-magnitudes. 


Filter 

E 

R,7(A) 

Aotab 

Am AC 

Am AD 

Amsc 

AmsD 

FWHM 

B 

5 

- 

0.19 ± 0.00!";“" 

-0.92 ± 0.00!°;“° 

-2.31 ±0.00!°;“! 

-1.11 ±0.00!°;“ 

-2.50 ± 0.00!“;“ 

0.30 

B 

7 

0.96 ± 0.04 

-0.01 ± 0.00!“;°“ 

-1.06 ±0.01!°;“° 

-2.36 ± 0.00!°;“ 

-1.05 ±0.01!°;“ 

-2.35 ± 0.00!“;“ 

0.30 

B 

9 

1.00 ±0.04 

-0.05 ± 0.00!“;°“ 

-1.11 ±0.00!°;“° 

-2.56 ± 0.01!°;“[ 

-1.07 ±0.00!°;“ 

-2.52 ± 0.00!“;“ 

0.30 

V 

6 

- 

0.33 ± 0.00!“;“° 

-0.82 ± 0.00!°;“ 

-2.19 ±0.00!“;“* 

-1.15 ±0.01!°;“ 

-2.52 ±0.01!“;“ 

0.24 

V 

7 

0.90 ± 0.03 

0.09 ± 0.00!“;“° 

-1.01 ±0.00!°;“° 

-2.31 ±0.00!“;“ 

-1.10 ±0.00!°;“ 

-2.39 ± 0.00!°;°“ 

0.24 

V 

9 

0.98 ± 0.03 

0.04 ± 0.00!°;“° 

-1.07 ±0.00!°;“° 

-2.51 ±0.01!°;°* 

-1.11 ±0.00!“;“ 

-2.55 ±0.01!“;“ 

0.24 

R 

3 

- 

0.52 ± 0.00!“;“° 

-0.66 ± 0.00!°;“* 

-2.15 ±0.00!“;“ 

-1.18 ±0.00!“;“* 

-2.68 ±0.01!°;“ 

0.13 

R 

4 

1.00 ±0.02 

0.53 ±0.01!°;“° 

-0.66 ± 0.00!°;“* 

-2.15 ±0.00!“;“ 

-1.20 ±0.01!°;“* 

-2.68 ±0.01!“;“ 

0.13 

R 

6 

0.99 ± 0.04 

0.41 ± 0.00!“;“° 

-0.76 ± 0.00!°;“* 

-2.14 ±0.00!“;“ 

-1.17 ±0.00!°;“ 

-2.55 ± 0.00!°;“ 

0.13 

R 

7 

0.90 ± 0.04 

0.14 ±0.00!°;“° 

-0.96 ±0.01!°;“* 

-2.23 ± 0.00!“;“ 

-1.10 ±0.01!°;“ 

-2.38 ± 0.00!“;“ 

0.13 

R 

9 

0.95 ± 0.01 

0.10 ±0.00!°;“° 

-1.00 ±0.00!°;“* 

-2.43 ± 0.01!°;“ 

-1.10 ±0.01!°;“ 

-2.53 ±0.01!“;“ 

0.13 

1 

9 

- 

0.01 ± 0.00!°;“° 

-1.01 ±0.01!°;“* 

-2.29 ± 0.01!°;“ 

-1.02 ±0.01!°;“ 

-2.30 ±0.01!“;“ 

0.30 

J 

1 

- 

0.34 ± 0.00!°;“ 

-0.64 ± 0.01!°;“[ 

-2.17 ±0.02!°;“ 

-0.99 ± 0.01!°;“ 

-2.51 ±0.02!“;“ 

0.29 

y#5133 

8 

- 

-0.01 ±0.01!“;°“ 

-0.84 ± 0.01!°;“ 

-1.78 ±0.01!°;“ 

-0.84 ± 0.01!°;“ 

-1.78 ±0.01!“;“ 

0.21 

Fmw 

7 

- 

-0.28 ±0.01!“;“ 

-1.16 ±0.01!°;!’ 

-2.57 ± 0.05!°;“ 

-0.87 ± 0.01!°;J“ 

-2.28 ± 0.05!“;“ 

0.08 

tf#5209 

8 

- 

-0.21 ±0.01!“;°“ 

-0.76 ± 0.01!°;“ 

-1.94 ±0.05!°;“ 

-0.56 ± 0.01!°;“ 

-1.73 ±0.06!“;“ 

0.21 

Ks 

2 

- 

-0.39 ± 0.00!“;“ 

-1.00±0.01!°;“[ 

-2.46 ± 0.00!°;“ 

-0.61 ± 0.01!°;“ 

-2.06 ± 0.00!“;“ 

0.21 

k'#5337 

8 

- 

-0.41 ±0.01!“;“ 

-1.01 ±0.01!°;“* 

-2.28 ± 0.04!°;“ 

-0.60 ± 0.02!°;“ 

-1.87 ±0.04!“;“ 

0.21 



Fig. 1. Simulated point-like magnitude difference AmAB(thin 
solid line), AniAc (dotted line) and Amsc (bold solid line) as 
a function of the average fraction of the flux of the host galaxy 
included in the total point-like flux after deconvolution (host 
contamination). 


nated by the ring than the flux of the bright ones. Consequently, 
the flux ratios between brighter and weaker lensed images are 
underestimated with respect to the theoretical values. This ring 
contamination increases with the ring brightness (i.e. from the 
B band to the K band) and we call this the background effect. 


The second effect is related to the specific geometry of 
J1131 in the source plane. Indeed, according to the lens mod¬ 
els, the QSO core is located at only 0.016"(115 pc in the source 
plane; Paper I) from a macro caustic. Hence, the dusty torus and 
the most central part of the host galaxy (that contributes signifi¬ 
cantly to the NIR flux) cover a region in the source plane with a 
strong magnification gradient. Therefore, the amplifications of 
the lensed images will depend on the considered source radius 
even at very small, unresolved scales. Moreover, since the local 
magnification around image A is larger than around images B 
& C, the unresolved source radius will be smaller for image A 
and the amplification ratios will also be affected. We dub this 
effect the differential amplification effect. 

Although the latter effect depends on the radial distribution 
of the source light, which is itself a function of wavelength, we 
made simulations to quantify its trend. We first chose a circu¬ 
larly symmetric exponential light profile with a half-flux radius 
re - 1.9 kpc (= 0.25") to represent a typical quasar host galaxy 
(Sanchez et al. I2004> . We lensed that extended source, sam¬ 
pled on a grid of 0.00025"/pixel, with a Singular Isothermal 
Ellipsoid + shear model fitting our data (SIE-ty; Sect. 14.2> us¬ 
ing the gravlens software vl.06 (Keeton. |2001> . The image 
plane, sampled on a 5 mas grid, was convolved with a Moffat 
profile (EWHM=0.7") and rebinned to a scale of 0.16"/pixel to 
mimic the observations. This synthetic frame was then decon¬ 
volved using the MCS technique to a final seeing of FWHM 












Sluse D. et al.: Multi-wavelength and multi-epoch imaging of RXS J113155.4-123155. 


7 


= 0.32". We repeated this operation for several source radii in 
the range 0.005"< < 0.1 "and we computed the flux ra¬ 

tios Ia/Ib and Ia/Ic- The main result is that the source remains 
unresolved in the A, B and C images for r^rc < 0.015". The 
associated flux ratios are then nearly identical to the point-like 
source ones (Am ~ 0.05 mag). For larger source radii, the im¬ 
ages become resolved and these flux ratios decrease and con¬ 
verge towards 1. 

Finally, we simulated a more realistic source model with 
two components: an extended exponential profile (truncated to 
'"host = 0.1") lensed into a resolved background ring and a su¬ 
perimposed disk-like source (rgso = 0.01") lensed into unre¬ 
solved images. This enabled us to look for the joined influence 
of the background ejfect and the dijferential amplification effect 
on the observed magnitude differences between the unresolved 
images. We quantify the effect as a function of the host contam¬ 
ination (defined as the fraction of the total point-like flux com¬ 
ing from the host galaxy). The results are displayed in Fig.^ 

Because the QSO Spectral Energy Distribution (SED) is 
expected to be dimmer than the host SED in the NIR range, 
we may expect a stronger ring contamination in that range. An 
accurate estimate would require us to know the true SEDs, but 
we can already infer an approximate contamination from the 
existing data. Indeed, in the B band, the ring image is faint and 
we estimate the host contribution to be at maximum 5% of the 
total point-like flux. Using a mean color index B - K - 2.5 for 
the QSO and B-K - 5 for the host galaxy, we calculate that in 
the K band up to 35% of the total point-like flux is coming from 
the host. Reporting these values in Eig.^ and comparing them 
with the case of 0% contamination, we derive that the observed 
flux ratios should not deviate by more than 0.05 mag from the 
predictions for a point-like source in the B band. In the K band, 
A/wab and Am\c are reduced by several tenths of magnitude. 
On the contrary, AmBc does not vary significantly. 


4. Results 

4.1. Ring morphology 

The multi-wavelength images of J1131 confirm that the con¬ 
trast between the QSO and its host decreases with increasing 
wavelength, as expected for the host galaxy of a quasar located 
at z ~ 0.66. Indeed, the SED of QSOs is minimum at ~ 1-2 
pm (rest-frame) while in this range the dominant emission is 
coming from the host (e.g. Dunlop et al. IT^ and marginally 
from the dust torus (hot dust emission increases up to ~ 2 pm-, 
e.g. Elvis et al. 1 19941 Nenkova et al. l2002t . 

While some structures in the ring are barely detected on the 
ground based Ks deconvolved frames (Eig.|2 right), the HST 
NIC2 deconvolved image (Eig. |2] left) spectacularly unveils 
many details in the ring. Additionally, a putative companion 
object (X) near the lensing galaxy is more clearly detected. Its 
nature is investigated in Paper II based on ACS and NICMOS 
data and it is shown that the expected influence of X on the lens 
modeling can only be very weak. 


4.2. Lens modeling 

In this section we construct a fiducial lens model aimed at dis¬ 
cussing the observed flux ratios and estimating the time delays 
for J1131. This adopted lens model is a Singular Isothermal 
Ellipsoid (SIE) -i- external shear (y), which is generically used 
in the study of the flux ratio anomalies (Dobler & Keeton l2005l 
and references therein). It has also the advantage to be fully 
constrained by the astrometry of J1131 (i.e. the positions of A, 
B, C, D and G). Additionally, it takes into account the main 
characteristics of the deflector, namely the lens galaxy elliptic- 
ity and the external shear associated with the lens environment. 
We used the gravlens modeling package v 1.06 (Keeton, 
1200 I t constrained by the HST relative astrometry of the point¬ 
like images and of the lensing galaxy (Sect. 1.3. 2t . Since sys¬ 
tematic errors affecting the NICMOS data (error on the plate 
scale, on the distorsion, on the orientation and due to PSE un¬ 
dersampling) do not influence similarly all lensed images, we 
have included them in the position error bars. The model yields 
the following results : Be - 1.843", e - 0.162, 0e - -57.92°, 
y - 0.114 and 9y - -83.01°; where 6e is the Einstein angular 
radius, e is the ellipticity and 6e is its PA, y is the tidal shear 
and 6j (counted E of N) points towards the mass producing the 
shear. We should note that the shear does not point towards any 
obvious massive perturber in the held. On the other hand e and 
6e agree well with the light profile ellipticity and PA. of the lens 
galaxy (see Paper II). The reduced (for 1 d.o.f.) associated 
with this model amounts to 136 which is statistically unaccept¬ 
able. The disagreement between the model predictions and the 
observations is mainly caused by the lensing galaxy position 
^gai “ expected 0.03" South from its observed position 
(TablellJ. 

Assuming that we can ignore X in the modeling, we 
opted for the SIE-i-y model as our fiducial lens model^. This 
model leads to predicted time delays between C (leading) 
and the other multiple images : tcb = 0.3 d, tca = 1-3 
d and tcd = 126.9 d. The predicted flux ratios are re¬ 
ported in Table |3 We have also fitted an ensemble of non 
parametric asymmetric models (searched by a Monte-Carlo 
method) with the PixeLens publicly available software (Saha 
& Williams I2004t . This enabled us to derive a distribution of 
Tcd (each tcd being associated with a different model) in the 
range 80-400 days while tca may not exceed 5 days. 

4.3. Flux variability 

4.3.1. Temporal flux ratio variations 

Hux variations of the lensed images of J1131 have been ob¬ 
tained in Sect. EH Eigure E represents the corresponding 
lightcurves in the R band* *. The comparison of these ones with 

^ We also emphasize that SIE+y left us with only 1 d.o.f. More 
complex models thus need additional constraints (see Paper II). 

* One might suspect that the difference of response between the 
Special R filter (epochs 3, 4, 6, 7) and the Bessel R filter (epoch 9) 
biases the observed lightcurve. But the color correction only amounts 
to 0.02 mag for the spectrum of J1131 and should not significantly 
modify our results. 
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Fig. 3. Lightcurves of images B, A, C and D in the R band for 
epochs 3, 4, 6 , 7, 9. The x axis refers to the observing date 
in modified Julian day (i.e. Julian day-2400000). The y axis 
refers to the magnitude difference with respect to epoch 3 (i.e. 
OT, - OT 3 ). The dotted line Joining the points is plotted for the 
legibility of the graph. 


the lighcurves observed in the V and B bands (i.e. epochs 6 - 
7-9 in V and epochs 7 and 9 in B) indicates an agreement be¬ 
tween the lightcurves within 0.05 mag and often around 0.02 
mag. Lightcurves are also nearly insensitive to the choice of 
the regularization parameter (except in the R band where the 
maximum shift of the D lightcurve can reach 0.05 mag). 

Figure l^reveals a dimming of B and C by about 0.3 mag 
in an eleven month period. A dimming of D is observed with a 
possible delay of about 5 months. During this period, the flux 
of A changed by only ~ 0.1 mag. 


Because the time-delay between A-B-C is of the order of 
a few days, intrinsic flux variations of the source would pro¬ 
duce nearly simultaneous flux variations between these images. 
Therefore, the different behaviour of image A with respect to 
B and C suggests that a microlensing event occurs in image A 
(i.e. microlensing partially compensates the intrinsic flux vari¬ 
ations; scenario SI). Alternatively, B and C might be dimmed 
by two independent microlensing events producing nearly si¬ 
multaneously the same lightcurve for both images (scenario 
S2). For both SI and S2, the observed microlensing time-scale 
is compatible with the crossing time of a microcaustic by the 
source (e.g. Treyer & Wambsganss l2004V Finally, we note that 
the dimming of D five months (i.e. ~ tcd) after B and C sup¬ 
ports SI, but microlensing in D cannot either be excluded. 

4.3.2. Chromatic variations 

Three phenomena can a priori introduce chromatic effects in 
the observed flux ratios: (i) differential extinction produced by 
dust in the lensing galaxy and which is expected to mainly 
affect the blue bands; (ii) microlensing, which is stronger 
for smaller source sizes (e.g.Wambsganss & Paczvhski ll991l 
Wyithe et al. I2000t and therefore for the shorter wavelengths 
emitted in the most central parts of the accretion disk (in the 
NIR, the microlensing signal would also be swamped by the 
constant flux contribution from the host); (Hi) systematic er¬ 
rors due to the differential amplification and to the background 
effect which decrease by several tenths the modulus of the mag¬ 
nitude difference between two images from B to K fSect. ITT^ 
and Fig.Q. 

Figure 0] shows the magnitude differences between the 
point-like images as seen through different filters from the B 
to the K band, at epochs 8 (NIR range) and 9 (optical range). 
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B V R I J H K 


Fig. 4. Magnitude differences between A (or B) and the other 
point-like images in B, V, R, I as observed at epoch 9 and in 
J, H, K at epoch 8 (cf Table 1^. The dotted line, when visible, 
represents the flux ratio predicted by the SIE-i-y model, namely 
Aotab = -0.54, A/ttAC = -0.65, Amec = -0.11 and Am ad = -3.23. 
Icr error bars are displayed. 

Chromatic variations are low (i.e. < 0.15 mag) in the optical 
while variations up to 0.5 mag are observed within the NIR 
range. This first statement does not argue in favor of strong 
differential extinction between lensed images since this one 
should be stronger in the optical range (where effect Hi is negli¬ 
gible) than in the NIR. Differential reddening can thus be safely 
neglected in the following. 

We now interpret the chromaticity of Ahiab, Anise and 
Aitiac (magnitude differences involving D are not discussed 
because of the large unknown time delay for this image) be¬ 
tween optical and NIR wavelengths. Note that in the following, 
we discuss the chromatic variations for optical and NIR range 
separately when microlensing might have produced flux ratio 
variations between epochs 8 -NIR- and 9 -optical-; 

- AmAB: Figure El shows that Aiiiab ~ 0 mag in the optical 
and decreases down to -0.45 mag in the K band. Both mi¬ 
crolensing scenarii described in Sect. EO can reproduce 
the fact that | AmAB I is smaller in the J band than in the K 
band; either A is affected (SI) and is then de-ampUfied or 
B is affected (S2) and is then amplified. Note that we do 
not expect chromatic variation caused by phenomenon Hi 
since Aiiiab ~ 0 mag in the J band. Whatever the adopted 
scenario, the value of Ib/Ia observed in the K band is the 
least affected by microlensing and thus constitutes an upper 
limit on the true flux ratio Ib/Ia (Tabled. 

- AmBc- lAmBcl ~ 1 mag in the optical range and decreases 
to lAmBcl ~ 0.55 mag in the H and K bands. Under SI, this 
fading of Ib/Ic from the B to the K band can be explained 
by effect Hi (since |AmBcl differs significantly from 0 and 
B and C are not microlensed). Note that since tcb is very 
short and since the ring contamination is negligible in the 
B band, AmBc in the B band should correctly represent the 
macro-lens flux ratio under S1. However, this value (AmBc 
= -1.05 mag) is abnormally low compared to the value pre¬ 
dicted by the model (Am"“‘^''^ = -0.11 mag). Under scenario 


S2, several combinations of microlensing for B and C and 
effect Hi are possible. 

- AiriAc ■ lAmAcI ~ 1-05 mag in the optical range. It amounts 
to 0.84 mag in the J band, 0.76 mag in the H band and rises 
up to 1.01 mag in the K band. First under SI, we can ex¬ 
plain the observed decrease of |AmAcl between the J and H 
bands by the dominant ring contamination, while the sub¬ 
sequent increase of IAtoacI in the K band would be due to 
the drop of the microlensing de-amplification of image A. 
Under S2, the same interpretation remains valid except that 
the increase of |AotacI in the K band implies that image C 
is more amplified by microlensing at bluer wavelengths. In 
the latter case, its value in the K band should be a lower 
limit of the macro-lens flux ratio. This reveals an abnor¬ 
mally large value of |AotacI (compared to the fiducial model 
prediction) under scenario S2. 

In summary, our data first suggest the absence of strong 
differential extinction between the lensed images of J1131. 
Second, the observed chromatic variations of the flux ratios can 
qualitatively be explained by the combination of effects ii and 
Hi. More interestingly, this analysis also reveals that either A 
is de-amplified hy a microlensing event (SI) or B and C are 
both amplified by two independent microlensing events (S2). 
For both microlensing scenarii, we have explained in which fil¬ 
ter the flux ratio measurement should be the closest to the val¬ 
ues predicted by the “macro model”. As one can see in Tabled 
these “best” estimates of the flux ratios do not reconcile with 
the macro model predictions. We investigate how anomalous 
are the flux ratios in the next section. 

4.3.3. Anomalous flux ratios 

Anomalous flux ratios in lensed systems have caught little 
attention until the work of Mao & Schneider iIT 9^ who 
first realised that substructures could significantly modify the 
flux ratios in B1422-H231. Various authors tried to compare 
the amount of Dark Matter Halo substructures predicted by 
ACDM hydrodynamical simulations to the amount of sub¬ 
structures necessary in lensing models to reproduce observed 
flux ratios (e.g. Metcalf & Madau l2001l Chiba 120021 Dalai & 
Kochanek EOO^ Metcalf & Zhao lMHl Mao et al. lMHl Bradac 
et al. |2002| |2004| Kochanek & Dalai IMmI Metcalf |2005i. 
However, it is still unclear whether the amount of substruc¬ 
tures necessary in lensing models agrees or not with the hy¬ 
drodynamical simulation predictions (see e.g. KGP05 or Rozo 
et al. l2005l for a recent extensive summary of the situation). 

Indeed, from a lensing perspective, flux ratios are not only 
sensitive to substructures but also to m=3 and m=4 multipoles 
of the lens model (Fvans & Witt 120031 Moller et al. 120031 
Quadri et al. 120031 Kawano et al. I2004t . to the lens environ¬ 
ment (Oguri|2003, to the source size (Dobler & Keeton l2005l 
Inoue & Chiba 120051 Chiba et al. I2005t and to isolated dark 
matter clumps on the same line of sight (Chen et al. 120031 
Metcalf l2005al I2005bl Wambsganss et al. 120051 . From a hy¬ 
drodynamical simulation perspective, the expected amount of 
dark matter within the Finstein radius of a lens is not well 
known. Indeed, both the abundance of halo substructures (e.g. 





















































10 


Sluse D. et al.: Multi-wavelength and multi-epoch imaging of RXS J113155.4-123155. 


Table 5. Summary of the observed relative astrometry and photometry related to J1131 and comparison with model predictions. 
The observed relative positions with respect to image A (obs), obtained from NICMOS images, are those given in Table|2 The 
systematic errors on the relative positions have been quadratically added to the random errors affecting the measurements. Since 
the position of A is not fixed to the observed one but known with the uncertainty reported in Table|2 the modeled position of A is 
not (0,0). Flux ratios /,/Ij in the K band are reported when microlensing (scenario S1 or S2) affects image i or j and in the B band 
otherwise. Since microlensing is minimum in the K band, the flux ratio reported is a lower/upper limit on the “true flux ratio”. 
Flux ratios including image D are placed into brackets because they are not corrected for the unknown time delay. Additionally, 
a microlensing event may also occur for that image. 


ID (Parity) 

Ad cos(5(") 


A<5('') 


Ib/I; 


Ia/I; 



obs mod-obs 

obs 

mod-obs 

SI 

mod 

S2 

mod 

A(-) 

0.000 ± 0.003 

-0.002 

0.000 ± 0.003 

-0.002 

<0.68 

0.61 

1.00 

1.00 

B(+) 

0.032 ± 0.003 

-0.002 

1.184 ±0.003 

0.003 

1.00 

1.00 

> 1.46 

1.65 

C(+) 

-0.590 ± 0.003 

0.002 

-1.117 ±0.003 

0.003 

2.67 ±0.01 

1.11 

>2.53 

1.83 

D(-) 

-3.115 ±0.003 

0.006 

0.875 ± 0.004 

0.019 

[10.02 ±0.01] 

11.95 

[10.58 ±0.06] 

19.72 

G 

-2.027 ± 0.003 

-0.010 

0.607 ± 0.004 

-0.029 

- 

- 

- 

- 

X 

-1.936 ±0.011 

- 

1.100 ±0.011 

- 

- 

- 

- 

- 


Zentner & Bullock [20021 Sigurdson & Kamionkowski 
and their spatial distribution (e.g. Chen et al. Zentner & 
Bullock l2003l De Lucia et al. 120041 Oguri & Lee l2004l are un¬ 
certain. In particular the amount of dark matter clumps that sur¬ 
vives in the central part of a halo (and so within a lens Einstein 
radius) is still debated (e.g. Mao et al. 120041 Amara et al. 120041 
Trott & Melatos 120051 Maccio et al. 120051 . 

In any case, the first step of any analysis is to obtain se¬ 
cure flux ratios and to identify anomalies independently of the 
choice of a macro-lens model. For this purpose, the most reli¬ 
able method is based on the validity/violation of the so called 
“magnification sum rules”, originally introduced by Schneider 
& Weiss The first rule predicts that a pair of images 

created when the source crosses a fold caustic will have the 
same brightness. The second one is valid when an unresolved 
source lies near a cusp caustic and predicts that the quantity 
^cusp (Mao & Schneider ll998t : 

^ _ \Ba + Bb + Bc\ _ |1 - Ib/Ia - Ic/Ia\ 

Ia^aI + Ia<bI + Ia<cI 1 + ffi/fA + fc/^A 

converges towards 0 as the source gets closer to a cusp singu¬ 
larity and all the three images are merging. 

The latter relation appears to be more robust than the re¬ 
lation binding “fold images” (KGP03, KGP05). Cusp config¬ 
uration systems are thus particularly important for the study 
of flux ratio anomalies (KGP03, Bradac l2004t . Among those, 
only two systems possess radio flux ratios (B1422H-231 and 
B2045 -h 265'^). The two others (i.e. RXJ 0911-1-0551 and J1131) 
could only be observed in the optical range and may thus be 
contaminated by dust extinction and microlensing. Hereafter, 
we further discuss the case of J1131 for which Rcusp ~ 0.35 at 
the epoch of the discovery (e.g. KGP03). 

Although microlensing precludes the determination of the 
true flux ratios, we have shown in Sect. 14.3.'^ that reliable up¬ 
per/lower limits on these ones could be derived (Table|3. Using 
the flux ratios of Tabled we And that Rcusp ^ 0.036 whatever 

* We do not consider the hybrid fold/cusp case of B0742-I-472. 


the microlensing scenario. Although Rcusp might naturally devi¬ 
ate from 0 (KGP03), the observed value is fully consistent with 
no violation of the cusp relation. This does not automatically 
mean that effects of substructures are rejected. First, we only 
have a lower limit on Rcusp and second, “small” Rcusp values can 
be observed while flux ratios cannot be reproduced by realistic 
lens models (cf B1422H-231; KGP03, Bradac et al. E^ . 

As already noticed in Sect. 14.3.'^ Ib/Iq (IaIIq) does not 
agree with our fiducial model predictions under SI (S2). 
Remembering that substructures can only magnify a minimum 
of the arrival time surface and preferentially demagnify sad¬ 
dle points (e.g. Schechter & Wambsganss l2()()2l . we now look 
if substructures are able to explain that flux ratio. Under sce¬ 
nario SI, Ib/Iq - 2.67 while the predicted value is 1.11 
(Table |3. Both B and C are minima of the arrival time sur¬ 
face. Consequently, only a magnification of B by 0.95 mag is 
allowed by the substructure scenario (since C being a mini¬ 
mum cannot be de-amplifled by milli-lensing) to explain IbIIq- 
However, the value of /a /7c (7a /7c > 3.92) derived under SI 
would then imply that A should be magnified hy a factor >2.15 
by substructures to match the predicted value of Ia/Ic - 1-83 
and this is very unlikely because i- two substructures are re¬ 
quired, ii- the saddle-point image A must be amplified. 

Under S2, IaHc S 2.53 while the macro model predicts 
1.83. Since image C is located at a minimum of the arrival 
time surface and, thus, cannot be demagnifled by substruc¬ 
tures, the above flux ratio can only be explained if image 
A is amplified by substructures by more than 0.35 mag (i.e. 
2.51og(2.53/1.83)). This would also reproduce the observed 
value of 7a/7b. Nevertheless, since image A is a saddle-point, 
it has a low probability to be magnified by a substructure. 

Consequently, under SI, a single substructure in front of a 
lensed image cannot explain the observed value of 7b IIq. Under 
S2, a milli-lensing magnification of A might explain the ob¬ 
served flux ratios. However, magnification of a saddle point by 
a substructure is statistically quite unlikely. 
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5. Summary and conclusions 

We have presented direct optical and NIR imaging of the 
quadruply imaged quasar RXS J113155.4-123155 obtained 
with ground-based telescopes and with HST at various epochs. 
The HST NICMOS images unveil many details in the Einstein 
ring of this system as well as a putative companion object (X) 
lying at from the lensing galaxy. Accurate positions have 
been derived for the 4 lensed images, for X and for the lensing 
galaxy based on the HST N1C2 data and on the ground-based 
VLT observations. Both sets of measurements do agree within 
6 mas. A SlE-ty model has been fitted to the observed image 
configuration. Nevertheless, the positions of image D and of 
the lens galaxy do not match satisfactorily the observed ones. 
Additionally, we have used the SIE-ty model to derive the ex¬ 
pected point-like flux ratios when taking into account the ex¬ 
tended nature of the source. 

The images have been deconvolved using the MCS method 
known to preserve photometry. Hence, we derived for the first 
time relative and absolute photometry of J1131 from the opti¬ 
cal to the NIR range. Since the photometric measurements of 
the point-like images in J1131 are complicated by the super¬ 
imposed Einstein ring, we have thoroughly discussed the po¬ 
tential sources of photometric errors. We have shown that the 
MCS method allows to reach a photometric accuracy close to 
the photon noise limit in the B and V filters if a good quality 
PSE can be constructed. On the other hand, small systematic 
errors affect the flux ratios in the R band due to the Einstein 
ring. 

The small chromatic variations of the flux ratios between 
the B and R bands suggest a low level of differential extinction 
between the lensed images of J1131. Based on sparse photo¬ 
metric measurements obtained between November 2002 and 
April 2004, we have shown the likely presence of microlensing 
for image A together with intrinsic quasar flux variability up to 
0.3 mag. Microlensing of A is also supported by our Chandra 
observations of J1131 in the X-rays (see Paper III, Claeskens 
et al., in preparation; Blackburne et al. I2005t . Additionally, 
based on the chromatic variations of the flux ratios, we have 
shown that A should likely be de-amplified (scenario SI), in 
agreement with the fact that it is a saddle-point of the arrival¬ 
time surface'® (Schechter & Wambsganss l2002l Keeton l2003t . 
Alternatively, our data cannot rule out that two independent mi¬ 
crolensing amplifications (scenario S2) with roughly the same 
amplitude occur for both images B and C. 

We have used the Rcusp relation (Mao & Schneider, 119981 
Eq.(^) to investigate the evidence for substructures in the lens¬ 
ing galaxy. Once lower/upper bounds on the flux ratios are im¬ 
posed from our knowledge of microlensing occuring in this 
system, we find a lower limit on Rcusp close to zero. This shows 
that the small scale structure(s) identified in KGP03 from the 
Rcusp relation, was (were) indeed due to (a) microlens(es). 

During the referee process, Rozo et al. <20051 demonstrated 
that small scale structures, if they can be considered as linear per¬ 
turbations, might not preferentially demagnify saddle point images. 
However, large perturbations from point-mass microlenses are likely 
caused by non linear perturbers and the results of Rozo et al. <20051 
may not be valid in such a case. 


Nevertheless, the observed flux ratios do not agree with the 
model predictions. We argue that milli-lensing effects cannot 
provide the best explanation under S1 while under S2, obser¬ 
vations suggest that a massive substructure should magnify the 
saddle point image A. This is however a low probability sce¬ 
nario (Schechter & Wambsganss I20031 l. An alternative to the 
milli-lensing explanation is to introduce in the macro-lens po¬ 
tential higher order multipole terms (e.g. Evans & Witt lM)^ . 
This is described in Paper II, where we successfully added an 
octupole (m = 4) to the fiducial model and reproduced both 
the observed flux ratio Ib/Iq (under SI) and the observed rel¬ 
ative lens position. Therefore, we found that the conjunction 
of microlensing and of higher order terms in the lens potential 
could explain the flux ratios observed in J1131. The preference 
of a more complex macro model over the presence of substruc¬ 
ture is opposite to what Kochanek & Dalai <200411 have found 
for other lenses. Of course the solution of adding an octupole 
might not be unique and other models such as those includ¬ 
ing several sources of shear (e.g. Morgan et al. l2005> . or with 
two mass components (e.g. Dye & Warren l2005ll should be ex¬ 
plored. 
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